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Effects of the electron-electron interaction in the two- 
dimensional flux phase are investigated. We treat the half- 
filled Hubbard model with a magnetic flux 7r per plaquette by 
the quantum Monte Carlo method. When the interaction is 
small, an antiferromagnetic long-range does not exist and the 
charge fluctuation is different from that of the Mott insulator 
It suggests that the Mott transition occurs at finite strength 
of the interaction in the flux phase, which is in contrast to 
the standard Hubbard model. 
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I. INTRODUCTION 

There has been a proposal that an order takes place on 
a link in several interacting lattice-fermion systems. Es- 
pecially when the link order has a phase factor, it brings 
an effective magnetic field. Sometime the order can be 
topological in the sense that the phase factor itself is not 
a well defined order parameter but the flux character- 
izes the phase. One of such systems is the flux phase 
which was proposed to describe the ground state proper- 
ties of several interacting lattice-fermion systems e.g. the 
Hubbard model, the t-J model, and their variants [0-0. 
Recently there has been a resurgence of interest in the 
flux phase and evidences are accumulating for its reality. 
For example, it was revealed a bond order takes place in 
the one-dimensional extended Hubbard model at half fill- 
ing Q which can be understood by the one-dimensional 
analog of the flux phase. Further, it has been discussed 
that a hidden topological order exists in cuprates, which 
shares some aspects with the flux phase M, In general, 
however, the flux phase competes with other instabilities 
e.g. superconductivity, antiferromagnetism, charge order 
and localization due to disorder (see refs. |I(]-[l3| for the 
effects of disorder on the flux phase). In particular, the 
flux phase often competes with superconductivity, which 
is a direct consequence from the SU(2) symmetry at half 
filling p|. We also note that the flux phase can emerge 
dynamically ^|J6) as well as in a static form [||. 

In this paper, we report effects of the interaction in 
the two-dimensional flux phase. We choose the two- 
dimensional Hubbard model with a magnetic flux <f> = tt 
per square plaquette and compare the results with the 
standard Hubbard model (0 = 0). In the standard Hub- 
bard model at half filling, it is believed that an infinites- 
imally small interaction drives the ground state to the 
Mott insulator, where a finite charge gap opens and an 



antiferromagnetic long-range order exists. This is con- 
sistent with the nesting argument in the weak coupling 
region. On the the hand, in the flux phase without the 
interaction, the density of state disappear linearly at the 
Fermi energy, which suggests that the structure of the 
low energy excitations is singular as compared with a 
simple Fermi liquid and the nesting instability is absent. 
Therefore one can expect an interaction-driven quantum 
phase transition from a singular quantum liquid (density 
of state is linearly vanishing without interaction) to a 
gapped insulator (Mott insulator). 



II. FLUX PHASE 

The flux phase is given by the ground state of the fol- 
lowing simple Hamiltonian, 
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where (J, k) denotes a nearest-neighbor link. The ampli- 
tude of tjk is constant but its phase factor tjk/\tjk\ — 
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e ie >" ^ satisfies a condition <j> = £ plaqucttc u jk 
a uniform magnetic flux per plaquette. The phase factor 
Ojk itself is not fixed but the flux <p is fixed, which is the 
gauge independent quantity. This Hamiltonian was pro- 
posed as an effective model (in the mean field level) of 
several correlated electron systems and discussed in many 
different contexts |l|-0,|ll] Oh One of the focuses was the 
stability of the flux state Jl4| |l6| . Following the discus- 
sion, the optimum, energy-minimizing, magnetic flux at 
half filling (this is the simplest case) is considered as 7r per 
square plaquette Jl4|,|l5|]. Furthermore, we note that Lieb 
gave some rigorous results for the stability in a general 
form 

At half filling, the low-lying excitations of the flux 
phase is described by massless Dirac fermions. There 
is a gauge freedom for the phase factor Ojk but let us fix 
them by choosing as tj +x j = (— 1)-M, tj +y j = t and oth- 
erwise zero where j = (j x , j y )&Z 2 , x — (1, 0), y = (0, 1). 
The energy bands in this gauge are given by 



EQs) = ±2tJcos 2 k x + cos 2 k y 



(2) 



w ±2t^J(k x - k x ) 2 + (k v - k y Y (i = 1, 2), (3) 

where (k x , k y )&[— n, it) x [0, n), k 1 = {k x ,k y ) — 
(7r/2,7r/2) and k 2 = (fc 2 ,fc 2 ) = (-tt/2, tt/2). Therefore 
the low-lying excitations are described by massless Dirac 
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fermions at these two gap-closing points and the density 
of states -D(e) near the Fermi energy vanishes linearly, 
D(e) oc |e|. The density of states -D(e) is singular and 
it leads to the suppression of the instability against the 
Mott insulator as discussed below. Note that the disper- 
sion is gauge dependent but the density of state is gauge 
independent. We focus on only the gauge independent 
quantity in this paper. 



III. MODEL AND METHOD 

We investigate effects of the interaction in the flux 
phase by the following Hamiltonian, 

U= ^j^jkCka + cl a t* k] c ja ) 

+ C/]T(n lT -l/2)(n a -l/2), (4) 

i 

where (j, k) denotes a nearest-neighbor link and U a on- 
site Coulomb repulsion. The geometry is set to be a two- 
dimensional square lattice and a periodic boundary con- 
dition is imposed. The grand-canonical ensemble is em- 
ployed and we put the system half-filled by the particle- 
hole symmetry. The \tjk\ is set to be constant (= t) and, 
based on the Lieb's theorem |lq ], the phase factor e l9jk is 
chosen so that the magnetic flux <fi is 7r per plaquette i.e. 
7r-flux phase (<j) = ^plaquette °jk =7i")- We always tries 
to compare the results of the flux phase (</> = ir) with 
those of the standard Hubbard model (0 = 0). It is to 
be noted that neither the translational nor time-reversal 
symmetry is broken in the 7r-flux phase. 

In order to study the system based on a non- 
perturbative approximation free method, the quantum 
Monte Carlo (QMC) technique is applied @Jl§]. We 
use the grand-canonical scheme at finite temperatures. 
Due to the particle-hole symmetry in the Hamiltonian 
(|J), the negative-sign problem does not occur. The sim- 
ulations were performed on a square lattice with a size up 
to N = 12 x 12 at a temperature down to T = 0.051 The 
Trotter decomposition is performed in the imaginary- 
time direction and the time slice is At ~ 0.10/i. We have 
checked that the systematic errors due to the Trotter de- 
composition are almost independent of temperatures and 
does not change the essential features after the extrap- 
olation. We have typically performed 500 Monte Carlo 
sweeps in order to reach a thermal equilibrium followed 
by 5000 measurement sweeps. The measurements are di- 
vided into 10 blocks and the statistical error is defined 
by the variance among the blocks. 

The Mott insulator is characterized by the following 
two features. One is a strong suppression of the charge 
fluctuation and the other is a presence of the strong anti- 
ferromagnetic spin correlation. In order to detect signals 
of the Mott transition, we have calculated the charge 
compressibility and the magnetic structure factor. The 
charge compressibility is defined by 
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where N e is the number of electrons and (3 an inverse 
temperature. The charge compressibility k measures the 
charge fluctuation directly. If the system has a finite 
charge gap, k shows a thermally-activated behavior and 
vanishes at zero temperature. The magnetic structure 
factor is given by 



S(q) 
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)((n iT -n a )(nj T -n^)). (6) 



If the system has an antiferromagnetic long-range order, 
S(tt, 7t) shows a diverging behaver as the temperature 
decreases. 



IV. RESULTS 

First let us discuss effects of the interaction on the 
charge compressibility. We compare the result with those 
of the standard Hubbard model to clarify the effects of 
the flux, that is, the structure of the low energy excita- 
tions. Fig|l] shows results of the charge compressibility 
k. Since the 7r-flux state at U/t = has gapless points 
in the Brillouin zone and the density of states D(e) near 
the Fermi energy vanishes linearly, D(e) oc |e|, the com- 
pressibility k does not show thermally-activated behav- 
ior. Even at U/t = 4, our data show that effects of the 
flux is still relevant and the charge gap is not well de- 
fined in the flux phase. This implies that the singular 
spectrum of the excitations seems to survive even with 
the interaction. Then this possible phase is clearly not 
a simple Fermi liquid but a singular quantum liquid. As 
the strength of the interaction increases, effects of the 
flux becomes irrelevant. If the interaction is sufficiently 
strong (U » t), the system becomes the Mott insulator 
with a finite charge gap which is the order of the inter- 
action. Therefore the results suggest the existence of a 
finite value of the interaction strength, U c , which sepa- 
rates from a gapless singular phase from the gapped one. 

Fig.|^ shows the antiferromagnetic structure factor 
S(tt, 7r) versus temperatures. For the standard half-filled 
Hubbard model, since the ground state has an antifer- 
romagnetic long-range order, S(tt, 7r) shows a diverging 
behavior as the temperature decreases (it saturates when 
the antiferromagnetic correlation length is longer then 
the lattice size). On the other hand, for the flux phase, 
the formation of the long-range antiferromagnetic order 
is not observed for U/t < 4. According to the spin-wave 
theory, S(tt, it) at the zero temperature increases with a 
lattice size as 



N 



m 

T 



+ 0(N- 1 / 2 ), 



(7) 



with m the staggered magnetization which is an order pa- 
rameter of an antiferromagnetic long-range order. Using 
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this relation, we try to obtain m 2 by plotting S(ir, ir) /N 
versus N -1 / 2 = L^ 1 . Fig.|| shows the plots for U/t = 4. 
The temperature is set to be T = 0.05i where the system 
reaches the zero temperature limit for the system size. 
For the standard Hubbard model, the data follow the re- 
lation (Q) and the extrapolation value is finite indicating 
the existence of the antiferromagnetic long-range order. 
On the other hand, for the flux phase, the relation (0) 
with m > does not hold i.e. there is no antiferromag- 
netic long-range order, which is in contrast to the stan- 
dard Hubbard model. When one discuss perturbatively, 
in the flux phase, the stoner instability is strongly sup- 
pressed due to the absence of the low energy excitations. 
The numerical results are consistent with this discussion 
at least in the weak coupling. 

Fig.^ shows the antiferromagnetic structure factor 
S(tt, 7r) for a variety of U ft. The antiferromagnetic corre- 
lation enhances as U/t increases. As noted above, when 
U/t is sufficiently large, effects of the flux become ir- 
relevant and one can expect that the antiferromagnetic 
long-range order appears. Due to the numerical diffi- 
culties, we can not perform simulations for stronger in- 
teraction (U/t > 10) regime. However, if the interac- 
tion is sufficiently strong (U^>t), the model (||) is es- 
sentially described by the antiferromagnetic Heisenberg 
model. Therefore the antiferromagnetic long-range order 
also may appear at some finite value of the interaction 
strength. 



V. DISCUSSION AND SUMMARY 

We have studied effects of the interaction in the flux 
phase. The Mott transition is focused using the quan- 
tum Monte Carlo method. Our results on the charge 
compressibility shows that effects of the flux is relevant 
for small U/t, while it becomes irrelevant when U/t is suf- 
ficiently large. The antiferromagnetism, which is charac- 
teristic of the Mott insulator, is also strongly suppressed 
in the weak coupling region. This is due to the structure 
of the low energy excitations in the flux phase, ft implies 
that the flux state with interaction leads to a new singu- 
lar phase for U < U c . This is in contrast to the standard 
two-dimensional Hubbard model. Effects of the doping 
is also an interesting future issue in connection with the 
competition with the superconductivity. 
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Fig. 1 




FIG. 1. Temperature dependence of the compressibility 
for a N — 10 x 10 lattice with interaction strengths U /t= 
(squares), 4 (circles), 6 (diamonds), and 8 (triangles). Solid 
symbols are for cj> — it and open symbols for <j> = 0. When the 
interaction is small (U/t < 4), the difference between cj> — 
and <f> = 7r is large over wide range of temperatures. 



Fig. 3 




FIG. 3. Extrapolation of antiferromagnetic long-range or- 
der. The dashed line is a least-squares fit to the data for 
4> = 0. For = 0, the points extrapolate to a finite value, 
indicating that the ground state has an antiferromagnetic 
long-range order. On the other hand, for 4> = it, S(tt, tv) versus 
1/L suggests the absent of an antiferromagnetic long-range 
order. 
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FIG. 2. The antiferromagnetic structure factor S(tv, n) as 
a function of temperature for (7/f = 4onaJV = 10x10 lattice. 
In the case of <f> — 0, S(tv, it) diverges at low temperatures due 
to the formation of the antiferromagnetic order. On the other 
hand, for <f> — n, S(tv, it) does not show a diverging behavior. 



Fig. 4 



.N=10x10 




<)>=7t j T , 






~\U/t = 8 T 




U/t = 6 


u ,. u/t = 4 


, , , , 1 , , , , 1 — 



°0 10 20 

P 



FIG. 4. The antiferromagnetic structure factor S(n, n) 
versus the inverse of temperature (5 for <j) — n on a N = 10 x 10 
lattice with interaction strengths U/t= 4 (circles), 6 (dia- 
monds), and 8 (triangles). The antiferromagnetic correlation 
enhances as U jt increase. 
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